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1 Introduction 



A new particle which can be interpreted as the Standard Model (SM) Higgs boson 
with the mass of around 126 GeV has been discovered at both the ATLAS and CMS 
experiments [1,2]. The discovery has a great impact on the supersymmetric models and 
supersymmetry (SUSY) breaking scenarios, since the 126 GeV Higgs boson requires large 
radiative corrections from top-stop loops within the minimal supersymmetric standard 
model (MSSM). The Higgs boson mass of 126 GeV is explained either by heavy scalar tops 
of O(10 — 1000) TeV [3] or by a large left-right mixing of the scalar tops [4] 1 . A simple 
and plausible possibility is that SUSY particles are as heavy as O(10 — 1000) TeV, which 
can ameliorate the Polonyi/moduli problem [6] and SUSY CP/flavor problems, though the 
setup is disfavored by the conventional naturalness argument. 

Besides the naturalness argument, light superparticles of O(100 — 1000) GeV are 
suggested by the muon anomalous magnetic moment (g — 2). The precise measurement of 
the muon g — 2 at the Brookhaven E821 experiment [7] shows a deviation from the SM 
prediction after the latest improvements on evaluations of the hadronic vacuum polarization 
contributions [8], 

5a^ = a% xp - = (2 6 .i ± 8 . ) x HP 10 . (1) 

This corresponds to more than three sigma discrepancy, where the hadronic light-by-light 
contribution referred to Ref. [9]. Although there are discussions on evaluations of the 
hadronic light-by-light contribution (see e.g., [10] for a recent review), they are expected to 
be settled in future, e.g., by lattice calculations [11]. A large uncertainty in the hadronic 
vacuum polarization contributions would be reduced, particularly if a disagreement among 
the experimental data for the 27r-channel region could be resolved [8, 12]. Moreover, 
new experiments are prepared to improve the precision of the experimental value at 
Fermilab [13] and J-PARC [14]. Consequently, although the discrepancy (1) is not a 
definitive evidence of physics beyond the SM, it is interesting to explore new physics 
models in light of the muon g — 2 as well as the Higgs boson mass of 126 GeV. In the 
MSSM, the above deviation of the muon g — 2 suggests sleptons, charginos and neutralinos 
at O(100) GeV and large tan/3 [15]. 

For light superparticles of 0(100-1000) GeV, gauge mediated SUSY breaking (GMSB) 
models are attractive from phenomenological points of view [16]. They are free from the 
SUSY flavor problem, and dangerous CP violating phases do not appear in scalar and 
gaugino masses. However, the models have difficulty in explaining the Higgs boson mass 
of 126 GeV as long as superparticle masses are O(100 — 1000) GeV, because the left-right 
mixing of the stops is too small to enhance the Higgs boson mass. As a result, GMSB 



1 Models with a large left-right mixing of the stops are constrained by the experimental bound from the 
branching ratio of the b — > S7 decay. See, e.g., Ref. [5] for the case of mSUGRA models. 
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scenario requires an extension to explain both the Higgs boson mass of 126 GeV and the 
muon g — 2. 

Recently, various extensions of the gauge mediation models have been studied for this 
purpose: the vectorlike-matter-assisted GMSB (V-GMSB) [17,18] 2 , additional U(l)' gauge 
symmetry imposed on the Higgs fields [5], a large left-right mixing of the stops induced 
by introducing a Higgs-messenger mixing [21], and splitting F-terms between colored and 
non-colored messengers [22] (see also [23]). Among the models, we will study the V-GMSB 
model in this paper. The model can enhance the Higgs boson mass easily by extra top-like 
matter multiplets in a similar manner to the top-stop contributions [24-26]. It has been 
shown that the model can explain both the Higgs boson mass of 126 GeV and the muon 
g — 2 simultaneously in a vast range of parameter space [17, 18] 3 . 

In this paper, we explore LHC status of V-GMSB in light of the latest SUSY search 
results of the ATLAS and CMS experiments. SUSY signatures of the V-GMSB model at 
the LHC have been briefly discussed in Refs. [18,20]. In this paper, the latest LHC results 
are applied to constrain the parameter space of V-GMSB, and we will show whether the 
region where the muon g — 2 is explained is still valid or not. Since the SUSY contributions 
to the muon g — 2 are suppressed as the superparticles are heavier, the V-GMSB model 
has an upper bound on the soft mass scale, e.g., on the gluino and squark masses. Thus, 
the model can be tested by direct searches for the superparticles. We will also improve 
evaluation of the superparticle mass spectrum. Relevant corrections from the vectorlike 
matters will be taken into account. It will be shown that the parameter regions of the 
gluino mass of < 1000— 1100 GeV are excluded in the case when the neutralino is the next- 
lightest-superparticle (NLSP), and more severe constraint is obtained for the stau NLSP. 
Consequently, it will be found that the whole regions where the muon g — 2 is explained 
at the la level are excluded by the LHC at 95% confidence level (CL). The region where 
the muon g — 2 is explained at 2a is still viable. 

This paper is organized as follows. In section 2, we will introduce the V-GMSB model, 
and discuss the Higgs boson mass, the mass spectrum, the muon g — 2, and constraints 
on the model. In section 3, analysis methods to study the collider searches for the 
superparticles will be described, and the latest LHC results will be applied. The main 
results are summarized in Fig. 2 and 3. Section 4 is devoted to conclusion and discussion. 

2 V-GMSB Model 

2.1 V-GMSB model and Higgs boson mass 

Let us explain the vectorlike-matter-assisted GMSB, or V-GMSB, briefly. The model 
with the GMSB boundary condition introduces a pair of the vectorlike matters in the 

2 The V-GMSB model was also studied in Refs. [19,20], where the muon g — 2 was not taken into account. 
3 Also, the model realizes the perturbative gauge coupling unification. 
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representations of 10 = (Q', U', E') and 10 = (Q 1 , U', E') of the SU(5) GUT gauge group. 
These extra matters generally have Yukawa couplings with the Higgs boson supermultiplets. 
The superpotential and soft SUSY breaking terms (below the SUSY-breaking scale) are 
respectively given by 

W = TUmssm + Y'Q'H U U' + Y"Q'H d U' + M Q ,Q'Q> + M W U'U' + M E ,E'E', (2) 

and 



r _ /"MSSM 
-t-soft — -'-soft 

+ m 2 Q ,\Q'\ 2 + m 2 Q ,\Q'\ 2 + m 2 uf \U'\ 2 + m 2 D/ \U'\ 2 + m 2 E ,\E'\ 2 + m 2 E ,\E'\ 2 

+ (A'Y'Q'HJJ' + A"Y"Q'H d U' + h.c.) 

+ (B Q ,M Q ,Q'Q' + Bjj.Mjj.U'U' + B E ,M E ,E'E' + h.c.) . 

The up- type (down- type) Higgs field is denoted as H u (Hd). We assume that the mixings 
between the vectorlike fields and MSSM fields are suppressed, though there must be small 
but finite mixings in order to avoid cosmological difficulties due to stable exotic particles. 
Such mixings induce decays of the extra matters, which will be briefly discussed in Sec. 3. 

The vectorlike top-like quarks, and t' 2 , and their superpartners give radiative 
correction via the Yukawa interaction Y'Q'H U U' to the Higgs boson mass in a similar 
manner to the top-stop loop corrections. This contribution is approximately expressed 
as [26] 



. 2 3Y'V 



ln M 2 eV 1 mjA 5 



M 2 \ A' 2 M 2 \ 1 A'^ 



+ T7o 1 



(4) 



.,,. s / , M 2 J M 2 \ 3M 2 J 12 Ml 

where the approximations of Mqi = Mjji{= M E ), tuq, = m 2 Q, = mfj, = m 2 j / (= Mj — M F ) 
and Bqi = B\ji = have been applied. In the numerical analysis, we include the full 
one-loop contributions. 

The radiative correction (4) becomes effective when the Yukawa coupling Y' is close 
to unity. Unless Y' is suppressed at a high scale, e.g., at the GUT scale, this is 
likely to be realized at low-energy scale, since it has a quasi-infrared fixed point at 
Y' ~ 1 [26]. On the contrary, the other Yukawa coupling Y" should be suppressed 4 , 
since the interaction Y" Q' HJLJ' decreases the Higgs boson mass [18]. If the couplings 
take the values of (Y',Y") ~ (1,0), the Higgs boson mass of 126 GeV is explained with 
Mqi , Mjji ~ 0(100 - 1000) GeV and M s ~ 0(1000) GeV. 

In the V-GMSB model, the soft SUSY-breaking terms are induced radiatively via 
messenger particles as in the ordinary GMSB scenario. The superpotential of the messenger 
sector is 

W mess = (M D + T D e 2 )^ D ^ D + {M L + T l 9 2 )VlVl, (5) 



4 The suppression of Y" can be explained in a model with the PQ symmetry [27], where Mq/ >( // jB / 
fi ~ 0(1) TeV is also related to the PQ-symmetry breaking scale. 
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where Mo (-Ml) is the super symmetric mass, and J~d(J~l) is the SUSY breaking F- 
term of the colored (non-colored) messengers. In the following, only a pair of the 5 + 5 
representation is included, which corresponds to the messenger number N mess = 1. The 
messenger scale is taken as M mcss = Md = M-l, and the SUSY breaking F-terms are 
as Td = Tl- Thus, the soft SUSY-breaking parameters are characterized by the two 
parameters, M. mcss and A = J^d/Md = Fl/Ml, at the messenger scale 5 . 

One of the characteristic features of the V-GMSB model is a size of the Higgsino mass 
parameter \i at the SUSY scale. It becomes very large compared to the usual GMSB 
models. This is because the vectorlike matters contribute to the renormalization group 
equation (RGE) of the up-type Higgs mass squared, which is approximately given by 



dm 2 jr dm 2 Tj 



d In Q d In Q 



+ i y ' 2 K' + <> + m k + Kl 2 )> ( 6 ) 

MSSM 



where Q is the renormalization scale. The complete formula at the two-loop level is found in 
Ref. [17]. The second term comes from the vectorlike matters, which gives a large negative 
contribution to the up-type Higgs mass squared rn 2 H ^ . It is as large as the stop contribution 
for Y' ~ 1. Hence, the electroweak symmetry breaking conditions require a large value 
of fj,. Consequently, the next-to-lightest neutralino x® and the lightest chargino xf tend 
to consist of the Wino, and the mass of the lighter stau is likely to be small, because the 
left-right mixing of the stau mass matrix, which is proportional to /utan/3, becomes large. 
Especially when tan j3 is large, the lighter stau becomes the NLSP. These features can be 
found in the mass spectrum discussed below. 



2.2 Mass spectrum 

In this subsection, we discuss the mass spectrum of the V-GMSB. In the numerical analysis 
we utilized SDFTSUSY 3 . 3 [28] and FeynHiggs 2 . 9 [29] to calculate the mass spectrum of 
the SUSY particles and the Higgs sector. They are modified to take the effects from the 
extra vectorlike matters into account. 

A typical mass spectrum is displayed in Fig. 1, where the V-GMSB result is compared 
to that of the ordinary GMSB model. Here and hereafter, we assume a common SUSY- 
invariant mass for the vectorlike fields, My = Mq> = Myi (= Mg/). As mentioned in the 
previous subsection, the heavier chargino xf and the heavier two neutralinos X3,4 mainly 
consist of the Higgsino, and are much heavier than the ordinary case. Also, it is found that 
the squarks become heavier. This is because the beta function of the strong gauge coupling 
constant is zero at the one-loop level, and the coupling stays large during the RGE. On 

5 The relation, Md = Ml and To = Tl, does not hold generally at the messenger scale due to the RG 
evolution, even if we take the universal condition, Mo = Ml and To = Tl, at the GUT scale. However, 
A is almost unchanged during the RG running, and hence, the soft mass parameters are not significantly 
changed. 
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Figure 1: The mass spectra of the GMSB model (left) and the V-GMSB model (right). The 
GMSB parameters are (A, M mess , tan (3, N mess ) = (140 TeV, 10 6 GeV, 20, 1) in both cases. 
The SUSY-invariant mass of vectorlike fields are set as Mqi = Mu' = Me> = 1 TeV for the 
V-GMSB model. The masses of vectorlike fermions (scalars) are labelled by r', b' , and 2 
(^1,25 &i,2> and *i,2,3,4)> respectively. 



the contrary, a ratio of the gaugino masses are less affected by the extra matters, since 
the ratios of Mi/ on, where Mi are the gaugino masses, are fixed during the RGE at the 
one-loop level. The masses of the scalar vectorlike matters are close to the top and bottom 
squarks up to the SUSY-invariant mass, because they have the same quantum numbers as 
corresponding squarks. The masses of the fermionic vectorlike matters are determined by 
the SUSY-invariant mass with the electroweak symmetry breaking corrections for t' . 

Let us explain how the mass spectrum is evaluated in SDFTSUSY and what kinds of 
modifications are applied. First of all, the program estimates the SM gauge couplings 
a a (Mz) and Yukawa couplings Yi(Mz) at the scale of the Z-boson mass Mz [30]. They 
are evolved upwards from Mz to the messenger scale M mess by solving the V-GMSB RGE, 
namely, including the contributions from the vectorlike matter supermultiplets. In the 
numerical calculations, the S0FTSUSY package is modified to include the two-loop RGE 
of the V-GMSB model, which can be found in [17]. It is also important to include 
threshold/self-energy corrections of vectorlike matters to the gauge coupling constants at 
the Mz scale, since they can give ~ 10% contributions especially to the strong coupling 
constant. They also affect the gaugino masses and the scalar masses squared generated at 
-^mess) which are respectively proportional to g 2 and g 4 . 

Next, the soft SUSY breaking parameters are provided by the messenger loops at 
M mess . The extra Yukawa couplings are also set at the scale, which are Y'(M mess ) = 1 
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and y"(M mcss ) = 0. The gauge, Yukawa and soft parameters are evolved by the V- 
GMSB RGE down to the SUSY scale, Mstjsy, which is determined by the stop masses as 
Msusy = V m *i m *Y 

The masses of the superparticles are evaluated including the whole one-loop corrections 
to the self-energies within the MSSM. The pole mass of the gluino also receives a correction 
from vectorlike matter loops. However, this turns out to be around a few GeV, since 
the masses of the vectorlike matters are close to Mstjsy, and their contribution to the 
self-energy are relatively small. Similarly, contributions of the vectorlike fields to the 
electroweak gaugino masses are safely neglected. 

Then, the Higgs potential is investigated, which determines the Higgs boson mass and 
the fi parameter. The MSSM part of the Higgs boson mass is evaluated at the two-loop 
level by the FeynHiggs package [29], while that from the vectorlike matters is estimated 
with the one-loop level effective potential [26], where all one- loop contributions from the 
vector like matters are taken into account. The two contributions to the Higgs boson mass 
are added in quadrature. The two-loop contribution of the vectorlike matters can shift 
the Higgs boson mass by ~ 1 — 10 GeV, study for which is reserved for future work. On 
the other hand, the /j, parameter is less affected by the vectorlike matters except for the 
RG evolution (6), because the tree-level contribution to fi is quite large in the electroweak 
symmetry breaking condition. 

There are several updates compared to the previous works in Ref. [17, 18], where the 
SuSpect2 [31] package was used to estimate the superparticle masses. In the old analyses 
the threshold/self-energy corrections to the gauge coupling constants from vectorlike 
matters were not included, which yield ~ 10% difference in the gauge couplings and the 
soft parameters for the GMSB scenario. Also, in Ref. [17,18], the fj, term contribution 
was discarded for simplicity in the estimation of the Higgs boson mass. In the present 
work, those contributions are included, which shift the Higgs boson mass by a few GeV. 
In addition, the self-energy correction from the MSSM particles to the stau mass is not 
fully included in SuSpect2. In the V-GMSB model, since the left-right mixing is relatively 
large, the lighter stau mass receives a sizable correction from the self-energy correction, 
leading to an O(10) GeV shift for a large tan/3. The mass spectra in the current analysis 
differs from those in Refs. [17,18] as well as Ref. [20]. 

2.3 Muon g-2 

In the V-GMSB model, the leading contributions to the muon g — 2 come from the smuon- 
neutralino and/or muon sneutrino-chargino loops as [15] 



Aa At (chargino) ~ 



sgn(//M 2 )tan/3, 



Aa At (neutralino) ~ 




sgn( / uMi) tan f3 + . . . , 



(7) 



24 ™soft 
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where m so f t is a typical scale of the soft parameters and the parameter. It is evident 
that the SUSY contributions are enhanced when tan /3 is large and m so ft is small. In 
the following, fx > is chosen, since the discrepancy (1) is positive. The vectorlike 
matters contribution to muon g — 2 is negligibly small. In the numerical analysis, the 
SUSY contributions to the muon g — 2 is evaluated by FeynHiggs, which includes two- loop 
contributions of the MSSM particles. 

2.4 Vacuum stability 

One of the most severe constraints on the V-GMSB model is the vacuum stability 
bound [32] . When tan (3 is large, the large trilinear coupling of the staus and the Higgs boson 
can generate charge breaking global minima and destabilize the electroweak symmetry 
breaking vacuum [33]. Given M mcss and A, tan/3 is tightly bounded from above by the 
condition that the lifetime of the electroweak symmetry breaking vacuum is longer than the 
age of the universe. In Ref. [32], a fitting formula given in Ref. [34] was utilized. Recently, 
however, the fitting formula has been revised as follows [34]: 

//tan/3 < 213.5^/m^mf R — 17.0(m^ 3 + mf R ) 

+ 4.52 x 10 -2 GeV -1 (m £s - m fR ) 2 - 1.30 x 10 4 GeV, 

where m^ 3 and m,f R are the soft masses for the left- and right-handed stau, respectively. 
Moreover, also very recently, the vacuum stability bound has been reanalyzed including 
effects of a radiatively-corrected tau Yukawa coupling [35]. The results in Ref. [35], however, 
cannot be directly applied to the V-GMSB model, because the masses and the mixing angle 
of the staus in our model are different from those in Ref. [35]. In addition, the thermal 
transition may change the bound [36]. In this paper we show the bound of Eq. (8), and 
also draw a bound which is weakened by 10% as a reference. 

2.5 Numerical Result 

Before moving on to the LHC constraints on the V-GMSB model, let us show the numerical 
results. The current status of the V-GMSB model is summarized in Fig. 2. (See also Fig. 3). 

In Fig. 2, contours of the Higgs boson mass and the muon g — 2 are drawn in the 
(rrig, tan /3) plane for the messenger scales, (a) M mess = 10 10 GeV, (b) 10 s GeV, and (c) 
10 6 GeV. Here, m g is the gluino pole mass, and tan/3 is a value at the scale, Q = Mz- 
In this paper, we use the strong coupling constant of a s (Mz) = 0.1184, and the top pole 
mass of m t = 173.5GeV. As discussed in Sec. 2.4, the vacuum stability bound is imposed 
by Eq. (8). In the figure, it is drawn by the blue dot-dashed line. The blue double-dotted 
long-dashed line corresponds to the vacuum stability condition which is weakened by 10%. 
The black solid and dashed lines show the LHC constraints, which we will discuss in detail 
in Sec. 3. In the gray shaded region, SOFTSUSY fails to calculate the stau mass. Note that 
such a parameter region is experimentally excluded, as long as the NLSP is long-lived. 
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In the green bands in Figs. 2 (a)-(c), the SM-like Higgs boson mass takes a value 
of 125 — 126 GeV for reference values of SUSY-invariant mass of the vectorlike matter, 
My = 1 TeV and 1.2 TeV. The Higgs boson mass increases as the vectorlike matter becomes 
lighter. The Higgs boson mass of 125 — 126 GeV can be realized in the whole parameter 
region of Fig. 2 by changing the vectorlike matter mass. For fixed My, the Higgs boson 
becomes lighter when the gluino mass is smaller, because the radiative corrections to the 
Higgs potential from (vectorlike-) top squarks decrease. On the other hand, if the messenger 
scale is higher, the (vectorlike-) top squark masses become larger due to the RG correction, 
and the Higgs boson mass of 125 — 126 GeV is realized for smaller gluino mass. 

Although the Higgs boson mass depends on My, masses of the MSSM superparticles 
and the SUSY contributions to the muon g — 2 are less sensitive to My. We have checked 
that the MSSM superparticle masses change only by < 2% when My is varied from 500 GeV 
to 1 TeV, and the muon g — 2 changes by less than 1%, correspondingly. In the figures, we 
fix My = 1 TeV in calculating all the quantities except for the Higgs boson mass. 

In the yellow (orange) region in Figs. 2 (a)-(c), the muon g — 2 discrepancy, Eq. (1), 
is explained at the 2a (la) level. As shown by Eq. (7), the region is sensitive to both the 
gluino mass, i.e., the soft mass scale, and tan/3. It is stressed that the soft mass scale 
has an upper bound in order to explain the muon g — 2 discrepancy. Also, too large tan (3 
spoils the vacuum stability condition. The bound becomes tighter as M mess increases. This 
is because the soft scalar masses are enhanced by RGE compared to the gaugino masses, 
which leads to suppression of the SUSY contributions to the muon g — 2 for a fixed value 
of the gluino mass. 

As we shall see in the next section, the LHC search for the SUSY particles depends on 
the species of the NLSP. In Figs. 2 (a)-(c), the NLSP is the neutralino (the stau) below 
(above) the light blue line, respectively. As M mess increases, the RG correction lifts up 
the scalar masses relative to the gaugino masses, and thus, the stau NLSP region becomes 
smaller. For the LHC search, the masses of the gluino, squarks, and the stau are particularly 
important. For illustration, we show in Fig. 2 (d) the masses of the lighter stau and the 
lightest squark among the first two generations by green and red lines, respectively, for 
M mcss = 10 6 GeV. The squark mass is almost independent of tan /3, since they are governed 
by the strong interaction. On the other hand, the stau mass, which determines the LHC 
constraint in the parameter region with the stau NLSP, depends both on the gluino mass, 
i.e., the soft mass scale, and tan/3 through the tau Yukawa coupling. 

3 LHC Limit 

Let us move on to the main topic of this paper. Constraints on the V-GMSB model from 
the LHC experiments will be explored. The current status of the superparticle searches at 
the LHC is applied to the V-GMSB model. 
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Figure 2: Contours of the Higgs boson mass, the muon g — 2 and the LHC constraints in 
the V-GMSB model are shown for (a) M mess = 10 10 GeV, (b) 10 8 GeV, and (c) 10 6 GeV. 
In the green bands, the SM-like Higgs boson mass is 125 — 126 GeV for the masses of the 
vectorlike matter, My = ITeV and 1.2 TeV. In the yellow (orange) region, the muon g — 2 
discrepancy is within 2a (lcr). Below the light blue lines, the NLSP is the neutralino, 
whereas the NLSP is the stau above them. The LHC constraints are shown by the black 
solid lines without theoretical uncertainties, where the regions left to the lines are excluded. 
The black dashed lines include theoretical uncertainties: 35% error of the production cross 
section is adopted for the neutralino NLSP, and 2% error in terms of the stau mass is 
for the stau NLSP. The regions above the blue dot-dashed lines are constrained by the 
vacuum stability condition, (8), and the bound which is weakened by 10% of (8) is drawn 
by the blue double-dotted long-dashed lines. The gray shaded regions are excluded by 
the failure of the SOFTSUSY calculation. In figure (d), masses of the stau and the lightest 
squark among the first two generations are shown respectively by green and red lines for 
M mess = 10 6 GeV. 
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In the model, the gravitino is the LSP, and collider signatures are determined by species 
and the lifetime of the NLSP. As discussed in the previous section, the NLSP is either the 
neutralino or the stau in the V-GMSB model, which are realized below and above the light 
blue lines in Fig. 2. In this paper, we concentrate on the case with a quasi-stable NLSP 6 . 
If a quasi-stable neutralino is the NLSP, searches for the events with no lepton, multi 
jets and missing energy give the most stringent constraint. On the other hand, long-lived 
slepton searches are the relevant search for a quasi-stable stau NLSP. These searches will 
be discussed in the following subsections. 

In the analysis, the mass spectrum is generated with customized SOFTSUSY 3 . 3 [28] as 
explained in the previous section, and passed to SUSY-HIT 1 . 3 [37] to calculate the decay 
pattern of the SUSY particles and the Higgs bosons. Prospino 2.1 [38-40] is utilized 
to calculate the next-to-leading order cross section. The CTEQ6L1 and CTEQ6 . 6M parton 
distribution functions (PDFs) [41] are used in the cross section calculation; the former is 
also used in Monte Carlo simulation for the event generation. 

Before proceeding to details of the LHC study, let us mention collider searches for the 
vectorlike quarks. In the V-GMSB model, the Higgs boson mass puts an upper bound 
on the masses of the vectorlike matters [18] 1 . Those particles are searched for in the 
LHC experiments. For example, the lighter top-like quark t[ has a lower mass bound of 
m? > 656 GeV if it decays exclusively into bW [42]. However, their decay pattern crucially 
depends on the mixing pattern between the vectorlike matters and the MSSM matters, 
may have other decay channels such as t'-y — > qZ and t' t — > qh, and then the constraints 
becomes much looser [18,42,43]. Therefore we will not discuss the searches for the vectorlike 
quarks anymore in the following. 

3.1 Neutralino NLSP 

For the region with a neutralino NLSP, we interpret the ATLAS result of the search for 
the superparticles in events with no lepton, 2-6 jets and missing energy, with data from 
yfl = 8TeV of f£ = 5.8 fb" 1 [44], as a constraint for the V-GMSB model. The CMS 
collaboration provides a similar bound [45], which is expected to result in similar constraints 
on the model. 

We employ Monte Carlo simulation to obtain constraints on the V-GMSB model. 

6 The NLSP could decay promptly if the gravitino is lighter than 0(1) keV. Then, collider signatures 
will change, such as events of di-photon with large missing energy for the neutralino NLSP case, and those 
of multi-taus with large missing energy for the case of the stau NLSP. Also, small violations of the R-parity 
can lead to a prompt decay of the NLSP. These cases require independent studies on LHC constraints, and 
are reserved for future works. See also discussion in Sec. 4. 

7 A precise evaluation of the upper bound on the vectorlike matter masses requires two-loop contributions 
of the vectorlike matters to the Higgs potential, which is beyond the scope of this paper. Note that the 
main results of this paper, the constraints from SUSY searches at the LHC, are almost independent of the 
masses of the vectorlike fields. 
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The kinetic distribution is calculated by MadGraph 5 [46] together with Pythia 6.4 [47] 
for simulating parton shower and initial/final state radiation, and Delphes 2.0 [48] for 
detector simulation. The events in pp — > gg,gq^*\q^q^ channels, which are relevant for 
the SUSY signals, are generated 8 . 

As for the object reconstruction, we follow the ATLAS procedure [44] up to the Delphes 
parameters. The Delphes original parameters for the ATLAS detector was used for the 
energy resolution and the calorimeter tower configuration. Jets are reconstructed using 
the anti-kt algorithm [49, 50] with AR = 0.4. Tracking efficiency for leptons (e and /i) 
is set 100%, which results in tighter lepton veto compared to the original analysis and 
thus conservative limits 9 . Only the jets with pp > 20GeV and |??| < 2.8, the electrons 
with pr > 20GeV and \r/\ < 2.47, and the muons with pr > lOGeV and \r/\ < 2.4 are 
considered. For overlap removal, jets are rejected if electrons are found within a distance 
of AR = 0.2, and then leptons within AR = 0.4 of any jets are discarded. Here, pr is a 
missing transverse momentum, AR = \J (Arj) 2 + (Acp) 2 is a distance between two objects 
in the (??, <f) coordinate, and r\ and <fi are pseudo-rapidity and azimuthal angle around the 
beam direction, respectively. 

Our analysis is employed with the same definition of the 12 inclusive signal regions 
(SRs) as the original paper [44], and we found that the "D-tight" SR among them yields 
the most stringent limit on the V-GMSB model (see Sec. 3.3). The selection procedure for 
the "D-tight" SR is defined as follows: 

1. having no lepton. 

2. having missing transverse energy $ T > 160 GeV. 

3. having at least five jets with pr > 60 GeV. In addition, the leading jet have pr > 
130 GeV. 

4. A0(jet, $i) > 0.4 for all of the leading three jets and > 0.2 for all jets with pp > 
40 GeV. 

5. $ T /rn§ > 0.15 and m^ c > 1700 GeV. 

Here, the effective mass m c ^ {m 1 ^) is defined as the scalar sum of Jvp and pr's of the 
leading five jets (all jets with pr > 40 GeV). 

The ATLAS search found no excess compared with the background estimation. The 
contribution from models beyond the SM (BSM) is constrained, which is shown as 95% CL 
upper limit on the number of the events as -/Vbsm < 6.0. This bound is interpreted as the 
constraint on the V-GMSB model. 

The production channels with i and 6 are neglected, whose production cross sections are less than a 
few % of the total cross section. 

9 Efficiency in lepton detection at the ATLAS experiment is estimated as ~ 90% and ~ 95% for electrons 
and muons, respectively, using y/s = 7TeV pp-collision data [51,52]. Including the lepton efficiency in the 
analysis results in an increase of the number of events which pass the following cuts; the increase is by 
~ 5%, which is smaller than theoretical uncertainties discussed later. 
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In order to validate the above procedure, we have checked that the ATLAS result for the 
mSUGRA model in Ref. [44] is reproduced in our simulation. The calculated distribution 
of the effective mass coincides with that of the ATLAS analysis. The exclusion limit on 
the mo — rrii/2 plane in our analysis is slightly looser than that of the ATLAS, but within 
la range of the theoretical uncertainty. 

3.2 Stau NLSP 

Searches for heavy long-lived charged particles is relevant for the case with a quasi-stable 
stau NLSP. We take the bound for the direct production of the lighter stau, pp — > f\f\, 
as the constraint on the V-GMSB scenario. The result is shown in the latest CMS report 
on the search for heavy long-lived charged particles with data at yfs = 7 TeV and fC = 
5.0 fb^" 1 [53]. It gives a lower bound on the mass of the lighter stau as 

m fl > 223 GeV (9) 

at 95% CL. Utilizing Prospino, we have checked that the cross section for the channel in 
the V-GMSB scenario is the same as that shown in the CMS report. 

The CMS collaboration also published the constraints on typical GMSB models, where 
the charginos and/or neutralinos are the main production channels [53]. However, it cannot 
be directly applied to the V-GMSB case, because the mass spectra are different between the 
two models (cf. Fig. 1). The difference results in a deviation of the velocity distribution of 
the NLSP stau. Once the reconstruction efficiency of the heavy long-lived charged particles 
is published, the production channels of the charginos and neutralinos can provide a more 
strict constraint. Otherwise, the inclusion of these channels will suffer from large systematic 
uncertainties. On the other hand, the bound for the direct production of the lighter stau is 
independent of the model, and it is conservative. The ATLAS collaboration also published 
a result on searches for the same signals [54] , but their search assumes direct productions of 
the whole sleptons, i.e., the production is not limited to the lighter stau. The reconstruction 
efficiency of the particle is not provided either. Since the mass spectrum of the V-GMSB 
model is different from the model analyzed in Ref. [54], we will rely on the CMS bound on 
the direct stau production in Eq. (9) in the following. 

3.3 Results 

Numerical results for the LHC constraints on V-GMSB are summarized in Fig. 2. Let 
us start from the neutralino NLSP region, i.e., below the light blue lines. The black 
solid lines represent the LHC constraint by searches for the jets with missing energy, and 
the black dashed lines show the limit with theoretical uncertainties, for which we adopt 
±35% errors in the production cross section. The uncertainties mainly come from errors 
of renormalization/factorization scales and a choice of the PDFs [55]. The left side of 
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(A,M mcss ,tan/3) 


(140TeV,10 6 GeV,20) 


(150TeV,10 10 GeV, 20) 


m-g 


1116 GeV 


1146 GeV 




1813 GeV 


2151 GeV 


production channel 


99 


M 


qq 


99 


gq 


qq 


cross section a (fb) 


5.54 


2.75 


0.238 


4.49 


0.702 


0.015 


acceptance A 


0.0720 


0.174 


0.189 


0.0858 


0.161 


0.138 


N = a x A x 5.81b- 1 


2.31 


2.77 


0.261 


2.23 


0.655 


0.0120 


N (total) 


5.3 


2.9 



Table 1: The cross section, the acceptance, and the expected number of events shown 
for the V-GMSB model points, (A, M mess , tan /3, M v ) = (140 TeV, 10 6 GeV, 20, 1 TeV) and 
(150 TeV, 10 10 GeV, 20, 1 TeV). The gluino mass, trig, and the right-hand squark mass of 
the first generation, m^ R , are also displayed. See Fig. 1 for the full mass spectrum of the 
former point. The acceptances are calculated for the "D-tight" SR, where 5000 events are 
generated at each model point. 

the black lines are excluded. Consequently, the gluino mass is required to be larger than 
about 1100 GeV for M mess = 10 6 GeV, if the theoretical uncertainties are not included. The 
theoretical uncertainties can shift the mass bound by < 100 GeV. When M mess is larger, the 
exclusion becomes weaker for the fixed gluino mass. This is because the squarks become 
heavier, and the production cross section of the squark, especially that of pp — > gif*\ 
becomes smaller. When the messenger scale is as large as 10 10 GeV, the bound becomes 
m g > 1030 GeV. 

The main production channels of SUSY events are pp — > gg, gq, qq. Although the 
squarks are relatively heavy compared to the gluino as shown in Fig. 1 and Fig. 2 (d) , they 
are not decoupled from the productions. Since the squark masses are almost independent 
of tan/3, so are the LHC exclusion lines in Figs. 2 (a)-(c). For illustration, in Table 1 
we show the cross section of each channel at two model points, (A, M mess ,tan/3,My) = 
(140 TeV, 10 6 GeV, 20, 1 TeV) and (150 TeV, 10 10 GeV, 20, 1 TeV). As can be seen in Fig. 2 
(a) and (c), they are close to the LHC exclusion limit. The acceptance and the expected 
number of events are also given in the "D-tight" SR. The acceptance is defined by a ratio 
of number of events that pass the cuts to the generated number of events. It is found that 
the channels of pp — > gg and gq comparably contribute to the SUSY searches for a low 
messenger scale (i.e., for a relatively light squark), while pp — > gg dominates for a high 
messenger scale (heavy squark). However, the contribution of the pp —> qq channel is small 
in the vicinity of the exclusion limits in Fig. 2 (a)-(c). 

The "D-tight" cut requires five hard jets and large missing energy. It is easy to 
understand the origin of the four hard jets: the gluino is lightest among the colored 
superparticles, and it generates at least two hard jets when it decays into lighter 
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Figure 3: Lower bound from the LHC on the gluino mass and the stau mass when the 
NLSP is the neutralino (left) and the stau (right), respectively. In the yellow (orange) 
region, the muon g — 2 discrepancy can be explained within 2a (lcr). 

superparticles. The origin of the fifth jet can be understood as follows. In the channel 
of pp —> gq, a hard fifth jet (or more jets) is generated by decays of the squarks into the 
gluino. In the pp — > gg channel, the largest fraction of events which pass the "D-tight" 
cut have additional hard jet(s) from initial/final state radiations (ISR/FSR) (c.f. [58]) 10 . 
The ISR can yield additional hard jet(s) in the pp — > gq channel as well. In the rest of the 
events, the fifth jet comes from decays of and r which are generated by the top quark 
decay and X\{X2) tt ^ ttxi fro m the colored superparticle decay, respectively. 

We have also checked that, in the mSUGRA model analyzed in Ref. [44], the parameter 
region which has squark and gluino masses similar to the V-GMSB model is mainly excluded 
by the "D-tight" SR. This result may support our conclusion that the V-GMSB model with 
a neutralino NLSP is mainly constrained by the "D-tight" SR among the 12 SRs of Ref. [44]. 

In the stau NLSP region, i.e., above the light blue lines in Fig. 2, the black 
solid line indicates = 223 GeV, corresponding to the 95% exclusion by searches 

for the heavy long-lived staus via the direct production [53]. The CMS collaboration 
evaluated the theoretical uncertainties of 3 — 7% in the cross section calculations due to 
renormalization/factorization scales, a s and PDFs. In the figure, we assigned theoretical 
error of 8% for the cross section, which corresponds to 2% error for the stau mass 
bound, shown by black dashed lines. The left-side regions of the black lines are excluded. 

10 The hardness of the ISR and FSR jets might not be simulated correctly by Pythia. In order to resolve 
this uncertainty, we analyzed "D-tight" by means of the MLM-matching scheme [56] implemented into 
MadGraph with a shower k t clustering and avoiding double counting between the gluino and squarks [57]. 
We have checked that the Pythia calculation agrees with the MLM-matching result. 
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Consequently, the region of the gluino mass lighter than 1200 GeV is excluded, and the 
bound becomes more severe for larger tan {3. Since the theoretical uncertainty is small, this 
bound does not change so much even if it is taken into account. 

It is interesting to analyze the LHC constraint as a function of the messenger scale in 
each category of the NLSP. The results are shown in Fig. 3. The left and right panels 
correspond to the cases of the neutralino NLSP and the stau NLSP, respectively. For a 
given M mess , the muon g — 2 discrepancy can be explained at the la (2a) level when the 
gluino/stau mass is within the orange (yellow) region. The upper ends of these regions 
represent the upper bound on the gluino/stau mass in order to explain the muon g — 2 
at the la (2a) level. In the case of the neutralino NLSP (left panel), the upper bound 
is determined by the requirement that the neutralino is lighter than the stau. In fact, 
for a fixed value of the muon g — 2, the gluino mass is maximized when tan/3 is set on 
the light blue lines in Fig. 2. When M mess is as large as 10 10 GeV, the vacuum stability 
bound can give a more severe bound. However, the result is almost unchanged, as noticed 
from Fig. 2 (a). This upper bound should be compared with lower bound on the gluino 
mass from the LHC constraint, drawn by the black solid line. Here, the exclusion limit 
is taken at tan/3 = 20 as a reference, since the bound does not depend much on tan/3. 
The black dashed lines include the theoretical uncertainties. (See discussion above.) It is 
found that all the regions where the muon g — 2 discrepancy is explained at the la level 
are already excluded by the direct searches for the superparticles at LHC. The region that 
can explain the muon g — 2 at 2a is still viable, and it is expected to be covered by LHC 
at ^fs = 13- 14TeV. 

For the stau NLSP case, the upper bounds on the stau mass are shown in the right 
panel of Fig. 3. The upper ends of the yellow and orange regions are obtained by the 
requirement that the stau is the NLSP. According to Fig. 2 (c) and (d), for a fixed muon 
g — 2 the stau mass is maximized when it is close to the neutralino mass. On the other 
hand, the black solid line is the lower bound on the stau mass from the LHC. It does not 
include the theoretical uncertainties; they are taken into account by the black dashed lines. 
Although the vacuum stability condition can give a tighter bound for M mcss ~ 10 9-10 GeV, 
such a parameter region is already constrained by the LHC (see Fig. 2 (a)). Consequently, 
it is found that all the parameter region that can explain the muon g — 2 at la is excluded 
by the searches for the heavy long-lived charged particles at LHC, and it is expected to 
become more severe, for instance by analyzing the data at yfs = 8TeV. 

4 Conclusion and Discussion 

The discovery of the Higgs-like boson had a significant impact on the scenarios beyond the 
Standard Model, particularly on the SUSY models. If the muon g — 2 is a real signal of the 
SUSY, the V-GMSB model is one of the only few viable SUSY models after the discovery 
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of the 126 GeV Higgs-like boson. 

In this paper, the collider searches for the superparticles in the V-GMSB model are 
explored in light of the current LHC results. The V-GMSB model predicts upper bounds 
on the gluino mass, i.e., the superparticle mass scale, to explain the muon g — 2 discrepancy 
with realizing the Higgs boson mass of 126 GeV. The LHC signatures are characterized 
by the NLSP. In the V-GMSB, the candidates are the neutralino and the stau. It was 
shown that the parameter regions of the gluino mass of < 1000 — 1100 GeV are excluded 
in the case when the neutralino is the NLSP, and the constraint is more severe for the 
stau NLSP. In both cases, it was shown that the whole parameter regions where the muon 
g — 2 is saturated at the la level are already excluded by the searches for the jets with 
missing transverse energy or those for the heavy long-lived charged particles. When the 
neutralino is the NLSP, all the region favored by the muon g — 2 at 2a will be covered by 
LHC after the upgrade for ^/s = 13 — 14TeV, while the whole region of the muon g — 2 at 
2a is expected to be covered earlier for the stau NLSP case. 

As shown in this paper, it is rather difficult to explain the muon g — 2 at la for the 
long-lived NLSP cases. One of the remaining possibilities is that the stau is the NLSP with 
a decay length less than (9 (cm). This short decay length is realized when the R-parity is 
slightly violated or the gravitino is as light as 0(10—1000) eV. GUT breaking effects on the 
SUSY breaking (invariant) masses of the doublet and triplet messengers can also help to 
avoid the constraints from the LHC results and to enhance the SUSY contributions to the 
muon g — 2 by generating more hierarchical masses of the colored and non-colored SUSY 
particles. 

Let us finally comment on the cosmology of the V-GMSB model. Assuming R-parity 
conservation, the gravitino LSP is completely stable, and the NLSP becomes long-lived 
unless the gravitino is ultra-light. The late-time decay of the NLSP may spoil the success 
of the big-bang nucleosynthesis (BBN), and hence there are constraints on the lifetime 
and the abundance of the NLSP, which translate into upper bounds on the gravitino 
mass, depending on the species and the mass of the NLSP. In Fig. 2, the NLSP mass is 
bounded from above in the region where the discrepancy in the muon g — 2 is reduced 
to 2a level; m-NLSP < 270 GeV for both of the neutralino and stau NLSP cases, for 
-^mess > 10 6 GeV. The BBN constraints on the gravitino mass are correspondingly 
^ 0.4 GeV for the neutralino NLSP and m 3/2 < 4 GeV for the stau NLSP [59]. 
The gravitino LSP with a mass of 0(1) GeV can be the dominant component of the dark 
matter if the reheating temperature is as high as T R = 10 6 - 10 7 GeV [60] 11 . Then, 
the non-thermal leptogenesis [62] may marginally work. We should also note that the 
notorious inflaton- induced gravitino problem [63] , which excludes many inflation models in 

llr The BBN constraints and the gravitino abundance depend on the mass spectrum of the MSSM and 
vectorlike particles. In particular, since the SU(3) gauge coupling remains strong up to the high energy scale 
in V-GMSB model, the perturbation in the calculation of the gravitino production becomes troublesome [61], 
which requires an improved analysis to reliably calculate the gravitino abundance. 
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the gravity-mediated SUSY breaking scenario with 0(1) TeV gravitino, can be avoided in 
(V-)GMSB models. Furthermore, the Polonyi problem can also be avoided in the GMSB 
models 12 ; rather, the decay of the coherent oscillation of the SUSY breaking field (Polonyi 
field) may lead to an alternative viable cosmological scenario [65]. More detailed study of 
cosmological scenarios in the V-GMSB model is beyond the scope of the present work, and 
will be given elsewhere. 
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